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Abstract 

We diseuss the preseription for the Dirae matrix in dimensional regularization used in most 
seeond- and third-order QCD ealeulations of collider cross sections. We provide an alternative 
implementation of this approach that avoids the use of an explicit form of y^ and of its (anti-) 
commutation relations in the most important case of no more than one y^ in each fermion trace. 
This treatment is checked by computing the third-order corrections to the structure functions F 2 and 
gj in charged-current deep-inelastic scattering with axial-vector couplings to the Vk-bosons. We 
derive the so far unknown third-order helicity-difference splitting function APns ’ that contributes 
to the next-to-next-to-leading order (NNLO) evolution of the polarized valence quark distribution 
of the nucleon. This function is negligible at momentum fractions x > 0.3 but relevant at x 1. 


Dimensional regularization [fTU^ . i.e., the analytie eontinuation of the theory to a non-integer num¬ 
ber D of spaee-time ‘dimensions’ (see also Ref. Q for an introduetion), is the standard framework 
for higher-order ealeulations in gauge field theories ineluding Quantum Chromodynamies (QCD). 
For some semi-leptonie benchmark observables, e.g., in inclusive deep-inelastic scattering (DIS) 
and semi-inclusive e^e~ annihilation (SIA), the use of dimensional regularization requires pre- 
scriptions for dealing with the genuinely four-dimensional objects £^vpo> the totally antisymmetric 
invariant tensor in four dimensions, and the Dirac matrix = iYoYiY 2 y 3 = */4! 

The tensor enters in the projection of the respective hadronic tensors onto the structure 
functions F 3 and gj in DIS and the fragmentation function Fa in SIA, e.g., 

= ... - ^ FAx, e"), (1) 

where the x is the Bjorken scaling variable, x = j{2pq) with = —q^, and where we have 
suppressed all non-F 3 parts of W^v h also occurs in the helicity-difference projection of incoming 
gluons in partonic polarized DIS. The matrix enters via the axial-vector coupling of the W and 
Z bosons to the quarks as well as by the corresponding helicity-difference projection for quarks. 

In particular the issue of y^ has attracted a considerable amount of attention. The ‘canonical’ 
approach is that of Ref. [|2l in which the Dirac algebra, and hence the loop momenta, are split in 4- 
and {D — 4)-dimensional sets with 

{Y5.Y4 = 0 . A' = 0, 1,2,3, 

[Ys’Ya'] = 0 , otherwise, ( 2 ) 

where {a,b} and [a,b] denote the standard anti-commutator and commutator, respectively. While 
Eq. Q leads to a consistent procedure [|4||, it has some drawbacks: the occurrence of additional 
scalar products of {D — 4)-dimensional loop momenta and an intermediate violation of the axial 
Ward identity. This situation has triggered quite a few of alternative suggestions which we are 
unable to address in this brief note; the reader is referred to [|5]-j9l and references therein. 

Our focus will be on the scheme developed, on the basis of Ref. ifTOll . in Refs. IHTUT^ which 
is closely related to that of Refs. OIHl but avoids complicating the loop integrals. Consequently 
this scheme has been employed in almost all higher-order (next-to-next-to-leading order, NNLO, 
or next-to-next-to-next-to-leading order, N^LO) diagram calculations of splitting and coefficient 
functions in DIS [[TTl - l2n and SIA [|^[23]l . as well as for the determination of the NNLO QCD 
corrections to the cross section for the production of a pseudoscalar Higgs boson Il24] - l26ll . 

We have been lead to consider this issue by our work on polarized charged-current DIS, in par¬ 
ticular the generalization of some of Ref. iflTll to the third order, which facilitates the determination 
of the so far unknown NNLO splitting function APns , the longitudinally polarized analogue of 
in Ref. [UTll . In order to study more cases with more than one Yg at the three-loop level, we 
have redone the calculations of F 2 and Fi of Refs. [[T5l[28ll^ and of gj of Ref. [[T^ with axial- 
vector instead of vector couplings to the gauge bosons. In particular for the latter case it was useful 
to employ an algorithm which is equivalent to, but more efficient than, that of Refs. lfTTI - fT3]| . This 
alternative implementation may be useful for future higher-order calculations in QCD. 
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In most of the higher-order ealeulations mentioned above, the preseription of Refs. [[TTI - fT3]| . 
sometimes briefly referred to as the Larin seheme, has been implemented in the form 

iy^^y5)L = , (3) 


i.e., what is eontinued is the axial-veetor matrix, written with a speeifie order of the two matriees. 
Alternatively one ean use (as, e.g., in Ref. II251I 1 

Y5,l = (4) 


Both substitutions, e.g., via Eq. ([T]), lead to produets of two £-tensors whieh ean be evaluated in 
terms of the D-dimensional metrie tensor 5^ as 






8 ^ S'" 8 ^ S'" 

8^ XV sv sv 
^p 

6S 8P SJ 
SS 5- 8- 8- 


(5) 


The need to use the D-dimensional metrie on the right-hand side has been elearly established, at 
least for the type of ealeulations we are eonsidering here, in Ref. lfT3ll . 

This implies that the traees 


Tr (yv Jv 2 ... Yv 2 ,„-iYpY 5 ) , 


( 6 ) 


evaluated using Eqs. (|3]) and dH) are not identieal at m > 3, as the additional terms generated by 
Eq. dll) eaneel only at D = 4 due to the Sehouten identity. 


p(4) 5V1 I p(4) 5^1 ip(4) X^l I p(4) X^l I p(4) 5V1 _ 

‘^V3V4V5V6 ^V2 ' ‘^V4V5V6V2 ' ‘^V5V6V2V3 'JV4 ' ‘-V6V2V3V4 ' t-v2V3V4V5 '^Ve ~ 


(7) 


However, if the above asymmetrie non-Hermitian form of the axial-veetor matrix is replaeed by its 
symmetrie Hermitian eounterpart (as done in Ref. [|25]I 1. 

YpY5 2(Y^Y5 - YsYa-) , (8) 

then Eq. dH leads to exaetly the same results as Eq. dl]) for the traee db]). The situation is eompletely 
analogous if the preseriptions dll) and dll) are applied, for any m and n, to 


Tr (yvJ v 2 ■ • • Yv„JpJ 5 YpiY p 2 ■ ■ • Yp^Y^Ys) : (9) 

Eq. dl]) and Eq. dH) with dl]) lead to the same results. Cases with more Y 5 will be addressed below. 

Obviously the ineonsistent use of Eq. dll) without Eq. dS]) leads to wrong results in diagram 
ealeulations only in suffieiently eomplieated oases. Eor example, re-oaloulating the third-order 
oorreotions for [lT^l20l in this manner leads to the same results as Eq. dl]) for eaoh individual 
diagram inoluding its dependenee on the gauge parameter. On the other hand, wrong (and unfao- 
torizable, ef. Ref. lHH) results would be obtained for the polarized veotor-axialveetor interferenee 
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structure functions g (using the labeling conventions of Ref. Il^ l in which occurs not once, 
as for F 3 , but twice. 

While the calculation of the Dirac traces is not usually a limiting factor in higher-order calcula¬ 
tions, the introduction of additional matrices by Eq. ([3]) or Eq. ([4]) has sometimes been considered a 
drawback of the Earin scheme. Eor traces with one y^, the most important case in QCD calculations 
(e.g., the only one encountered in Refs. lfT3l - l^ l. this issue can be avoided by using algorithms 
which are completely equivalent and do not introduce any additional intermediate matrices. 

A procedure equivalent to, but faster than, using Eq. ([3]) is provided by the following steps: 

1. Write the one-y^ traces in the form ® without changing the order of the y-matrices. This 
can be viewed as using the cyclicity of the trace, or as reading it from this point, cf. Ref. [I 6 l|. 

2. Evaluate Eq. ® using 

Tr (yviyv2 • ■ ■ yV2„,-iyuY5) = ~^^Sviy2 ■ ■ ••?V2„,-5V2m-4^V2m-3V2m-2V2m-l/' 

± permutations of Vi...V 2 ,„_i . ( 10 ) 

Incidentally, this main step can be programmed in Eorm [l3TI - [3^ . for an extensive docu¬ 
mentation see [l34ll . in a very compact manner for any number of traces with one y^, viz 

repeat; 

id,once,G(ml?,?a,mu?,five) = distrib_(-2,3,Gl,G2,?a)*G(mu,five); 
id G2(mul?,mu2?,mu3?)*G(mu4?,five) = e_(mul,...,mu4); 
endrepeat; 

. sort 
repeat; 

if ( count(Gl,1) ); 

id,once,Gl(?a) = g_(l,?a); 

Tracen,1; 
endif; 
endrepeat; 

3. Eor traces with more than one y^, use Eq. ([3]) for all but one (special care is needed for more 
than two y^, see below), then calculate the resulting one-y^ trace according to 1 . and 2 . above. 

A corresponding algorithm equivalent to Eq. (|4l) can be implemented by changing 1. and 2. above to 

1. Input all axial-vector matrices in the form ([ 8 ]), then proceed as under 1. above. 

2. Evaluate the resulting traces, in which now y^ has no special role, as 

Tr (yviyv2 ■ ■ • Yv2,„-iyuy5) = ~^^8viy2 ■ • ■'?V2m-5V2m-4^V2m-3V2m-2V2m-l,u 

± permutations of Vj. • • V 2 ,„_i , (11) 

for which the central two lines of the above EORM implementation are changed to the simpler 
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id,once,G(ml?,?a,five) = distrib_(-2,4,Gl,G2,?a); 
id G2(mul?,...,mu4?) = e_(mul,...,mu4); 


Eq. (fTTI) has certainly been used elsewhere before; however we have not seen a elear diseussion of 
the ‘implicit-Y 5 ’ relations (fTO and (fTTI) to the ‘explicit’ prescriptions dH) and dll) in the literature. 

We now have four equivalent manners to evaluate traees with and should briefly address 
their effleiency: eomputing Eq. dS) for m = 7, i.e., with 14 y-matriees besides y^, requires about 
1.2 and 38 seconds, respectively, using d3]) and dH) with the internal trace algorithms of Eorm, but 
0.3 and 1.2 seconds with the shown implementations of Eqs. (fTO and (fTTI) on a Xeon E5-2667v2 
with 3.30GHz, using one core. The eorresponding numbers for m = 8 are higher by faetors of 
about 20. This sealing is the same as for the non-y^ ease, which is however faster by almost a 
factor of 8 than our fastest y^ implementation (fT^ . The eorresponding execution times for Eq. dl]) 
with rn = n = 5 (12 y-matriees besides the two y^) are, in the same order, 4.5, 740, 1.3 and 55 
seconds; the two faster methods again take longer by about a faetor of 20 for m = n = 6 . 

We now move to the application of the above y^ scheme in higher-order calculations, focusing 
on the best known (in general and to us) ease of third-order DIS in massless perturbative QCD. 
This seheme shares the second drawback of the ‘t Hooft/Veltman seheme dH), the violation of the 
axial Ward identity. This issue is less serious here than it may be in higher-order ealeulations in the 
electroweak theory; it is addressed by ‘eorreeting’ the axial eurrent by the renormalization factors 
Z 5 and Za determined to the third order in the strong coupling constant tts in Ref. IfTTI . 


Z 5 = 1 —as4C^-|-a, 


-1 2 


, ( 12 ) 


22 c} 


107 I 2 ^ 

-g- Z,pL,A -f- g Z-pUj: 


+ a' 


Ci?( - 2|e + 96C31+c/=ca - I 6 OC 3 1+C,Cil - + 


.2^ / 5834 


2147 


-|- Cpfij- 


62 32 


356 , 32 


27 - — ^3 ) +CACpnf[ + ~ ^3 ) + 81 CpU 


81 


52 


81 


(13) 


for D = 4 — 2z. These factors are expressed in terms of the renormalized eoupling normalized as 
as = tts/( 47 t), and we have employed the first two eoeffieients of the P-function of QCD [l35l - [38]l . 

o _ 11 ^2 n _ 34 ^2 10 ^ /I 

Po — 3 “ 3 7 Pl ~ 3 ^A “ 3 3 • ( 14 ) 

to write Eq. dT^ in a slightly more compaet form. 

On top of, or instead of, the multiplieation with Z 5 Za before performing the mass faetorization, 
a non-trivial faetorization-seheme transformation is required in the polarized ease for arriving at the 
splitting and coefficient funetions in MS for the helieity-dependent ease, see Refs. lfT9ll^[39l - l4^ . 
At N^LO this transformation is not fully known yet: the pure-singlet quark eontribution is missing. 

A seeond yet innocuous effect of using Eqs. dll) - dS]) (or any equivalent algorithm) is that all 
traees, ineluding those of the a® Born eontributions, reeeive an additional dependence on D. This 
dependence is faetorized and then removed in the projection on the strueture functions. E.g., the 
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well-known D-dependence in the projeetion on the strueture funetion F 3 , 



1 

^ (D-3)(Z)-2) 


^ 5 


p-q 


( 15 ) 


originates in the basie trace of with four other y-matrices and ~ (D — 2) (D — 3). This 

factor is analogous to the {D — 2)^^ in the Fj projection that arises from ypy^yP = {2 —D) y^. 

As mentioned above, Eq. ([3]) has been extensively used in higher-order QCD correction in cases 
where only one y^ occurs. On the other hand, we are not aware of a corresponding NNLO or N^LO 
calculation involving two occurrences y^ in either the same or different traces. The former case is 
more interesting and challenging; a good first example is a third-order calculation of the structure 
functions F 2 and F^, for typical forward-Compton diagrams see Fig.[Tl with an axial-vector instead 
of the vector coupling [l28l|29l|43l|44l to the gauge boson. 




Figure 1: Typical third-order Feynman diagrams for the two flavour classes contributing to quark- 
initiated charged-current DIS. Depending on the structure function, the boson lines at the top are 
replaced by or a combination of and p^Py, and the quark lines at the bottom by jp or y^y^ 
(in Schoonship notation). The vertices with p and V represent vector or axial-vector couplings. 


Since we are not looking for a new splitting or coefficient function in this calculation, it is 
sufficient to keep the full dependence on the Mellin variable N only to two loops and determine 
the third-order corrections for a few even-integer values of N. At this level the computation is 
straightforward and virtually automatic; we can use our old diagrams databases and employ our 
calculation and analysis programs with minor modifications. 

The projection on the axial-vector structure functions F 2 and Fi involves a different prefactor. 
The Bom-level trace (p/q are the quark/gauge-boson momenta, = 0) is now 


Tr(ypy^y 5 yp+,Y''Y 5 ) ~ {D-l){D-2){D-3){D-6)pq , (16) 

i.e., the projections involve an extra factor [(D— 1)(D —3)(6 —D)]~^ Taking this into account, and 
multiplying the results by (Z 5 Z^i)^ as given by Eq. (fT^ and (fTSI) before factorization, we obtain the 
same splitting functions and quark and gluon coefficient functions for F 2 and Fl as found before. 


(«)aa 

fq/g 


(x) 


(/i)vv 

fq/g 


for 


/ = 2, L at n = 1, 2, 3 . 


(17) 
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This demonstrates that there is no need to resort, as often done, to a fully anti-eommuting in 
traees with two y^ (whieh admittedly would lead to the right result here): the seheme eonsidered 
here ean be used also for these eases, at a usually tolerable eost in eomputing time. 

The ultimate ehallenge, in the framework of QCD eorreetions for structure functions at the 
lowest order in the electroweak theory, is provided by doing the same for structure function gj in 
polarized DIS. This calculation involves, in addition to that for F 2 , an £-tensor from the projection 
on gj, essentially the same as Eq. (fTSl) . and a y^ or £-tensor taking the quark or gluon helicity 
difference. The resulting contractions of four £-tensors have to be performed with special care as 
their order matters in the present case, unlike in four dimensions where the results can be shown 
to be the same by repeated application of the Schouten identity ([7]), cf. Ref. Il45l . 


In the case at hand, it is correct to pair the £-tensors from the axial-vector vertices (labelled 
and V in Fig.fl]). This is readily achieved in Form by using the built-in tensor e_ only for these, 
and to ‘protect’ the other two £-tensors by using a different notation until the other contractions 
and traces have been performed. The fastest implementation is to use Eq. dS]) for the axial-vector 
gauge-boson vertices together with Eq. (fTOl) . with e_ suitably renamed in the Form code shown 
below that equation. A three-fold application of Eq. ([3]) is also possible, if considerably slower. 
The symmetric implementations dll) and (fTTI) are yet less efficient; the four-fold application of dH) 
is prohibitively slow at the third order. Only now all four prescriptions consistently lead to 

Tr(y„T 5 Y;.T 5 Tp+,YvY 5 )e'’’'" ~ (D-2)(D-if(D-6)(pqf, (18) 


enabling us to verify, by diagram calculations. 


(n)aa , . _ (n)vv , . 

gi,q/g^ ^ 


at n 


1,2,3. 


(19) 


Finally we address the splitting functions APns for the polarized quark-antiquark differences 

Mr = a4-a4-(a4-a4) , (20) 

”/ 

A4 = E{a4-a 4} Pi) 

of helicity-dependent parton distributions. A/, = — f^~, where and f[~ represent the distri¬ 

butions of the parton i with positive and negative helicity, respectively, in a nucleon with positive 
helicity, and is the number of effectively massless flavours. For general reasons one expects also 
at NNFO, n = 2, a direct relation between the polarized and unpolarized non-singlet cases 

, ( 22 ) 

of which the right-hand side was calculated to NNFO in Ref. ifTTll . On the other hand, the difference 

API, = AF;, - AP„; (23) 

can only be determined by a diagram calculation. It is this calculation, via the two-y^ polarized 
vector-axial vector interference structure function gg (cf. Ref. [[30]l l that lead to our above consid¬ 
erations on yg. In particular, APjP^ is obtained from the flavour class //02 in Fig. [U where the 
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M^-bosons are not attached to the external quark line, for the helicity projection = JpJ^ and 

the structure function projection i.e., with the two entering in different traces. 

The resulting even-A^ Mellin-space expression reads 

APiP\N) = 16n^J'^^'^4fc,K(5_3(-20ri + 8ri2)+5i,_2(8ri-16ri2) 

+ 32ri S_2,i + 53 (6ri + 4ri2) + S_2 (8ri + 20ri^ + Sri^ - 4D^) 

+ Si (8ri - 14^^ - 42ri3 - -2D^)+D^^ (24) 


in the notation of Ref. ff 2 TI . i.e., with 0^ = {N + k)~^,r\= DqD^ and all harmonic sums [l46ll taken 
at argument N. The corresponding x-space result, in terms of harmonic polylogarithms [|47]I at 
argument x (also suppressed), is given by 

AP„?)^(x) = 16n^rf“^^'4i,,K((l-x)((24-20C2)Hi-8C2Ho,-i-2Hi,o,o 
— I6H0-1-1,0 + 8Ho,-i,o,o + 8Ho,o,-i,o) + (1 +Jc)( 16 H_yo — 52 ^ 2 H_i 
-8^2Ho,1 - 40 H_i _i,o + 36 H_i,o,o + 32 H_i,o,i + 12 Ho,0,0,1 - 4 Ho,1,0,0) 

- X (16 ^3 Ho + 8 Ho.o + 36 Ho,o,o - 8 Ho,o,o,o) - Hq( 1 + 24x) - Ho (6 - 74x )^2 

- Ho, 0(12 + 20 x )^2 + Ho,i (10 + 8 x) + Ho,-i,o (8 + 36x) + Ho,o,i (6 - 38x) 

-(10-8x)C2 + (6 + 88x)C3 + (25 + 15x)C4) , (25) 


which can be parametrized, with an accuracy of about 0.1% or better for 10 ^ <x < 0.95, by 

APiP\x) ^ nf{l-x) (-42.97 - 29.29 Lq + 179.1 + 117.8x - 385.5^2 + 15.94x^ 

+ xLo(8.818Lo + 460.8)+2.6811n(l-x))+0.0001ny5(l-x) , (26) 

where Lo = Inx and the artificial 5(1—x) contribution can be included to compensate the slightly 
lesser accuracy at very large x to improve the approximation for high-// moments and large-x 
convolutions. Eqs. (l24l) and (l25l) . together with our calculational verification of Eq. (l22l) from 
the even moments of at NNEO, cf. Ref. [|27]l . complete the determination of the third-order 
helicity-dependent splitting functions of which the main part was performed in Ref. Il2n . 

Eq. (I 25 ]) can be employed to determine also the odd moments, in particular 

APn?^''(A^= 1) = ^n^r^^dabcln, (23 - 12^2 - 16^3) ■ (27) 


Together with 


(A^=l) 

APi?”(iV=l) 


= Q(Q-2Q) (-134-12^2-8^3) , 


(28) 


= Cj{Cj,-2Cp) 
+ (Q — 2Cp^ 


145 

2 

1081 

36 


62^2 + 164 ^3 - 372 ^4 + 48 ^2^3 + 208 ^5 


, 245,, 3214,, , 1058,, ,„,,,, 

+ ^ ^2 --Q-^3 + ^4 - 48 C 2 C 3 - 112^5 


'^r''i/^26 44 y 448 y 68 ^ 

LFnf(LA-2Lp)[ 9 + 3 ^ 2 - 9 S3+ 3 ^4 


( 29 ) 
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- note the presenee of ^2 and the higher weight in the Riemann-^ funetion as eompared to the 
‘natural’ even moments, ef. Seetion 3 of Ref. Il49ll - this leads to the expansion 

APnP''(iV= 1) ^ - 0.00810a^ - (0.04075 -0.01850n^) + 0(a4) (30) 

in QCD, i.e., for = 3, Cp = 4/3 and d^^^^dabc/f^c — 5/18. For the normalization of 

the latter we had to ehoose between the earlier eonvention of Refs. ifTTlfT^fTTHTOll and that of 
Refs. [|20ll44]| and other more reeent artieles following Eq. (187) of Ref. Il48]l . We have done the 
latter; however the reader should be aware that the unpolarized eounterparts of Eqs. (l24l) - (l25l) 
were presented using the eonvention d^’^^^dabc/^c = 40/9 in Ref. ifTTl . 

Returning to Eq. (l30l) . we note that for Hj- = 4 Eq. (ITTI) provides almost two thirds of the 
eorreetion, whieh is aetually larger than the tiny part at normal seales; without it the eoeffieient 
of would only amount to 0.00061. At large-v is negligible though: it is suppressed 
by two powers of ( 1 —x) with respeet to AP^, with the leading large-x term the same as for its 
unpolarized eounterpart in Eq. (4.11) of Ref. IfTTlI with 32n;y(i“^‘^Ja^c/^c(2C2 — 3)(1— v)ln(l—x). 

0.4 


0.2 


0 


- 0.2 


-0.4 

0 0.2 0.4 0.6 0.8 1 lo"'’ 1010'' 

X X 

Eigure 2: The NNEO splitting funetions and for the polarized quark distributions 

(I 20 I) and (I 2 TI) . together with the previously unknown leading eontribution (l25l) to their differenee 
(I 23 ]) for three flavours, divided by 2000 ~ (47t)^ to eompensate for our small expansion parameter 
as = (Xs/ ( 47 t). Also shown, on the right, is the unpolarized eounterpart llETlI of Eq. (1251) . 

The situation is totally different at small x, as shown in Eig. [2) despite an only quadratieally 
logarithmie (negative) small-x enhaneement, 

APi?\x) = -l 6 nfd^^^dabc/n,{ 6 ^ 2 Ll + (l+ 6 ^ 2 )Eo + 0(1)} , (31) 

its eoeffieients are sueh that it overwhelms at x > 10 ^^ the (positive) small-x behaviour of AP^^ 
whieh ineludes terms up to In^x that are, due to Eq. (I22l) . given by Eq. (4.15) of Ref. [[T71 . 




10 '^ 10 '^ 10 
X 


8 















To summarize, we have diseussed some subtleties of the (multiple) use of the preseription of 
Refs. llTTlfT^ with a D-dimensionally eontraeted £-tensor [[T3]l in higher-order QCD ealeulations, 
and provided a proeedure that is eonsiderably faster than the algorithm mostly used so far and 
henee may be useful in some future three- and four-loop ealeulations in QCD. We have applied our 
findings to re-derive some third-order results in polarized and unpolarized deep-inelastie seattering, 

(2)s 

and to ealeulate the hitherto unknown NNLO splitting funetion AP^s (x) whieh eontributes to the 
evolution of the polarized valence quark distribution, thus completing the determination of the 
NNLO splitting functions for helicity-dependent parton distributions of hadrons. 

A Form procedure of our alternative implementation of the scheme of Refs. lfTTI - [r3]l . as well 

(2)s 

as Form and Fortran files of our results for APns be obtained by downloading the source 
of this article from http: //arxiv. org/ or from the authors upon request. 
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